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Abstract. Thermal photon production at mid-rapidity in Au-l-Au reactions at = 200 GeV is studied in the frame¬ 
work of a hydrodynamical model that describes efficiently the bulk identified hadron spectra at RHIC. The combined 
thermal plus NLO pQCD photon spectrum is in good agreement with the yields measured by the PHENIX experiment 
for all Au-l-Au centralities. Within our model, we demonstrate that the correlation of the thermal photon slopes with 
the charged hadron multiplicity in each centrality provides direct empirical information on the underlying degrees of 


freedom and on the form of the equation of state, s{T)/T^, 
the reaction. 

PACS. 12.38.Mh - 24.10.Nz - 25.75.-q - 25.75.Nq 


1 Introduction 

Numerical calculations of lattice QCD predict a transition from 
ordinary hadronic matter to a deconfined state of quarks and 
gluons when the temperature of the system is of the order of 
Tcrit ~ 0.17 GeV jT). The existence of such a phase transition 
manifests itself clearly in the QCD equation-of-state (EoS) on 
the lattice by a sharp jump of the (Stefan-Boltzmann) scaled 
energy density, e(7’)/7’^, at the critical temperature, reminis¬ 
cent of a first-order phase change*. The search for evidences 
of this deconfined plasma of quarks and gluons (QGP) is the 
main driving force behind the study of relativistic nuclear col¬ 
lisions at different experimental facilities in the last 20 years. 
Whereas several experimental results have been found consis¬ 
tent with the formation of the QGP both at CERN-SPS |2l and 
BNL-RHIC 13 energies, it is fair to acknowledge that there 
is no incontrovertible proof yet of bulk deconfinement in the 
present nucleus-nucleus data. In this paper, we present a de¬ 
tailed study of the only experimental signature, thermal pho¬ 
tons, that can likely provide direct information on the thermo¬ 
dynamical properties (and, thus, on the equation-of-state) of 
the underlying QCD matter produced in high-energy heavy-ion 
collisions. Electromagnetic radiation (real and virtual photons) 
emitted in the course of a heavy-ion reaction, has long 141.51 
been considered a privileged probe of the space-time evolution 
of the colliding system^, inasmuch as photons are not distorted 
by final-state interactions due to their weak interaction with 

* The order of the phase transition itself is not exactly known: the 
pure SU(3) gauge theory is first-order whereas introduction of 2-1-1 
flavours makes it of a fast cross-over type Q. 

^ Excellent reviews on photon production in relativistic nuclear col¬ 
lisions have been published recently I6t7i8l . 


of the strongly interacting matter produced in the course of 


the surrounding medium. Direct photons, defined as real pho¬ 
tons not originating from the decay of final hadrons, are emit¬ 
ted at various stages of the reaction with several contributing 
processes. Two generic mechanisms are usually considered: (i) 
prompt (pre-equilibrium or pQCD) photon emission from per¬ 
turbative parton-parton scatterings in the first tenths of fm/c of 
the collision process, (ii) subsequent y emission from the ther- 
malized partonic (QGP) and hadronic (hadron resonance gas, 
HRG) phases of the reaction. 

Experimentally, direct y have been indeed measured in 
PbH-Pb collisions at CERN-SPS = 17.3 GeV) 0. 

However, the relative contributions to the total spectrum of the 
pQCD, QGP and HRG components have not been determined 
conclusively. Different hydrodynamics calculations 
1121131141 require “non-conventional” conditions: high initial 
temperatures (7 q"“ > Tcrit), strong partonic and/or hadronic 
transverse velocity flows, or in-medium modifications of 
hadron masses, in order to reproduce the observed photon 
spectrum. However, no final conclusion can be drawn from 
these results due mainly to the uncertainties in the exact 
amount of radiation coming from primary parton-parton 
collisions. In a situation akin to that affecting the interpre¬ 
tation of high pt hadron data at SPS ca, the absence of 
a concurrent baseline experimental measurement of prompt 
photon production in pH-p collisions at the same ^/s and pj 
range as the nucleus-nucleus data, makes it difficult to have 
any reliable empirical estimate of the actual thermal y excess 
in the Pbi-Pb spectrum. In the theoretical side, the situation at 
SPS is not fully under control either: (i) next-to-leading-order 
(NLO) perturbative calculations are known to underpredict 
the experimental reference nucleon-nucleon y differential 
cross-sections below « 30 GeV fT^ (a substantial amount 
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of parton intrinsic transverse momentum kj CD, approx¬ 
imating the effects of parton Fermi motion and soft gluon 
radiation, is required mi (ii) the implementation of the extra 
nuclear kj broadening observed in the nuclear data (“Cronin 
enhancement” m resulting from multiple soft and semi-hard 
interactions of the colliding partons on their way in/out 
the traversed nucleus) is model-dependent 11201211221 and 
introduces an additional uncertainty to the computation of the 
yields, and (iii) hydrodynamical calculations usually assume 
initial conditions (longitudinal boost invariance, short ther- 
malization times, zero baryochemical potential) too idealistic 
for SPS energies. The situation at RHIC (and LHC) collider 
energies is undoubtedly far more advantageous. Firstly, the 
photon spectra for different centralities in Auh-Au Ga and in 
(baseline) p-tp 1241 collisions at ^ = 200 GeV are already 
experimentally available. Secondly, the p-tp baseline reference 
is well under control theoretically (NLO calculations do not 
require extra non-perturbative effects to reproduce the hard 
spectra at RHIC I24IB1 ). Thirdly, the amount of nuclear 
Cronin enhancement experimentally observed is very modest 
(high pt are barely enhanced in dn-Au collisions at -yiivw = 
200 GeV (26l), and one expects even less enhancement for y 
which, once produced, do not gain any extra kj in their way out 
through the nucleus. Last but not least, the produced system 
at midrapidity in heavy-ion reactions at RHIC top energies is 
much closer to the zero net baryon density and longitudinally 
boost-invariant conditions customarily presupposed in the 
determination of the parametrized photon rates and in the 
hydrodynamical implementations of the reaction evolution. 
In addition, the thermalization times usually assumed in the 
hydrodynamical models (t,/,„,„ < 1 fm/c) are, for the hrst time 
at RHIC, above the lower limit imposed by the transit time of 
the two colliding nuclei (Tq = 2R/YRi 0.15 fm/c for Auh-Au at 
200 GeV). As a matter of fact, it is for the hrst time at RHIC 
that hydrodynamics predictions agree quantitatively with most 
of the differential observables of bulk (“soft”) hadronic pro¬ 
duction below pr « 1.5 GeV/c in Auh-Au reactions IITTE^ . 

In this context, the purpose of this paper is three-fold. First 
of all, we present a relativistic Bjorken hydrodynamics model 
that reproduces well the identihed hadron spectra measured at 
all centralities in Auh-Au collisions at y/svw = 200 GeV (and, 
thus, the centrality dependence of the total charged hadron mul¬ 
tiplicity). Secondly, using such a model complemented with the 
most up-to-date parametrizations of the QGP and HRG pho¬ 
ton emission rates, we determine the expected thermal photon 
yields in Auh-Au reactions and compare them to the prompt 
photon yields computed in NLO perturbative QCD. The com¬ 
bined inclusive (thermaln-pQCD) photon spectrum is success¬ 
fully confronted to recent results from the PHENIX collabora¬ 
tion as well as to other available predictions. Thirdly, after dis¬ 
cussing in which pr range the thermal photon signal can be po¬ 
tentially identihed experimentally, we address the issue of how 
to have access to the thermodynamical properties (temperature, 
entropy density) of the radiating matter. We propose the corre¬ 
lation of two experimentally measurable quantities: the thermal 
photon slope and the multiplicity of charged hadrons produced 
in the reaction, as a direct method to determine the underlying 
degrees of freedom and the equation of state, s{T)/T^, of the 


dense and hot QCD medium produced in Auh-Au collisions at 
RHIC energies. 


2 Hydrodynamical model 

2.1 Implementation 

Hydrodynamical approaches of particle production in heavy- 
ion collisions assume local conservation of energy and momen¬ 
tum in the hot and dense strongly interacting matter produced 
in the course of the reaction and describe its evolution using the 
equations of motion of perfect (non-viscous) relativistic hydro¬ 
dynamics. These equations are nothing but the conservation of: 

(i) the energy-momentum tensor: = 0 with = (£ + 

p)u^u'^ — pg^'^ [where e, p, and m'' = (y,Yv) are resp. the en¬ 
ergy density, pressure, and collective flow 4-velocity fields, 
and g^''=diag(l,-l,-l,-l) the metric tensor], and 

(ii) the conserved currents in strong interactions: dfjjf = 0, 
with jf = HiU^ [where n, is the number density of the net 
baryon, electric charge, net strangeness, etc. currents]. 

These equations complemented with three input ingredi¬ 
ents: (i) the initial conditions (Eq at time Xo), (ii) the equation- 
of-state of the system, p(e,n, ), relating the local thermodynam¬ 
ical quantities, and (ii) the freeze-out conditions, describing the 
transition from the hydrodynamics regime to the free streaming 
final particles, are able to reproduce most of the bulk hadronic 
observables measured in heavy-ion reactions at RHIC II27I28I 
|29l. 

The particular hydrodynamics implementation used in this 
work is discussed in detail in IT 2 I We assume cylindrical 
symmetry in the transverse direction (r) and longitudinal (z) 
boost-invariant (Bjorken) expansion GqI which reduces the 
equations of motion to a one-dimensional problem but results 
in a loss of the dependence of the observables on longitudinal 
degrees of freedom. Our results, thus, are only relevant for 
particle production within a finite range around midrapidity^. 
The equation-of-state used here describes a first order phase 
transition from a QGP to a HRG at Tcrit = 165 MeV with latent 
heat^ Ae « 1.4 GeV/fm^, very similar to that used in other 
works frh . The QGP is modeled as an ideal gas of massless 
quarks {Nf = 2.5 flavours) and gluons with total degeneracy 
gQGP = few+ 7/8g,„„fa) = 42.25. The corresponding EoS, 
p = 1/38 —4/3B (B being the bag constant), has sound 
velocity cl = dp/de = 1 /3. The hadronic phase is modeled as 
a non-interacting gas of ~400 known hadrons and hadronic 
resonances with masses below 2.5 GeV/c^. The inclusion of 

^ The experimental 7t^ and dN/dy distributions at RHIC are 
Gaussians 1311 . as expected from perturbative QCD initial condi¬ 
tions 1321 . Thus, although there is no Bjorken rapidity plateau, the 
widths of the distributions are quite broad and within |y| < 2, devia¬ 
tions from boost invariance are not very large 1321 . 

^ Although the lattice results seem to indicate that the transition is 
of a fast cross-over type, the predicted change of Ae ~ 0.8 GeV/fm^ in 
a narrow temperature interval of AT ~ 20 MeV Q can be interpreted 
as the latent heat of the transition. 
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heavy hadrons leads to an equation of state significantly differ¬ 
ent than that of an ideal gas of massless pions; the velocity of 
sound in the HRG phase is ~ 0.15, resulting in a relatively 
soft hadronic EoS as suggested by lattice calculations l33l : 
and the effective number of degrees of freedom at Tc is g„f.a « 
12 (as given by = A5s /see later). Both phases 

are connected via the standard Gibbs’ condition of phase 
equilibrium, Pqgp{Tc) = PHRa{Tc), during the mixed phase. The 
external bag pressure, calculated to fulfill this condition at Tc, 
is B fv 0.38 GeV/fm^. The system of equations is solved with 
the MacCormack two-step (predictor-corrector) numerical 
scheme 03 with time and radius steps; 5f = 0.02 fm/c and 5r 
= 0.1 fm respectively. 

Statistical model analyses of particle production in nucleus- 
nucleus reactions ca provide a very good description of the 
measured particle ratios at RHIC assuming that all hadrons are 
emitted from a thermalized system reaching chemical equilib¬ 
rium at a temperature Tchem with baryonic, strange and isospin 
chemical potentials p,. In agreement with those observations, 
our specific hydrodynamical evolution reaches chemical freeze- 
out at Tchem = 150 MeV with ps = 25 MeV (as given by the 
latest statistical fits to hadron ratios IShl l. and has ps = Pi = 0. 
For temperatures above Tchem we conserve baryonic, strange 
and charge currents, but not particle numbers, while for temper¬ 
atures below Tchem we explicitly conserve particle numbers by 
introducing individual (temperature-dependent) chemical po¬ 
tentials for each hadron. The final differential hadron dN/dpr 
spectra are produced via a standard Cooper-Frye ansatz Ell 
at the kinetic freeze-out temperature (Tfi, = 120 MeV) when 
the hydrodynamical equations lose their validity, i.e. when the 
microscopic length (the hadrons mean free path) is no longer 
small compared to the size of the system. Unstable resonances 
are then allowed to decay with their appropriate branching ra¬ 
tios EH. Table I summarizes the most important parameters 
describing our hydrodynamic evolution. The only free parame¬ 
ters are the initial energy density Eq in the center of the reaction 
zone for head-on (impact parameter b = 0 fm) Auh-Au colli¬ 
sions at the starting time Xq, and the temperature at freeze-out 
time, Tfo- 


2.2 Initialization 

We distribute the initial energy density within the reaction vol¬ 
ume according to the geometrical Glauber^ prescription pro¬ 
posed by Kolb et al. Bol . Such an ansatz ascribes 75% of the 
initial entropy production in a given centrality bin, io (^) Uo soft 
processes (scaling with the transverse density of participant 
nucleons Npart{b)) and the remaining 25% to hard processes 
(scaling with the density of point-like collisions, Ncoii{b), pro¬ 
portional to the nuclear overlap function TAA{b)): 

s{b) = C ■ (0.25 • Npartib) + 0.75 • Nccii{b)), (1) 

^ The density of participant and colliding nucleons are obtained 
from the nuclear overlap function 7 aa(^) computed with a Glauber 
Monte Carlo code which parametrizes the Au nuclei with Woods- 
Saxon functions with radius R = 6.38 fm and diffusivity a = 

0.54 fm EH- 


where C is a normalization coefficient chosen so that we 
produce the correct particle multiplicity at b = 0 fm. For each 
impact parameter, we construct an azimuthally symmetric 
hydrodynamical source from the (azimuthally deformed) 
initial Glauber entropy distribution, by defining a coordinate 
origin in the middle point between the centers of the two 
colliding nuclei and averaging the entropy density over all 
azimuthal directions. We then transform ^o{b) 

This method provides a very good description of the measured 
centrality dependence of the final charged hadron rapidity 
densities dNch/dry measured at RHIC as can be seen in 
Figure [2 Note that in our implementation of this prescription, 
we explicitly added the contribution of the particle multiplicity 
coming from hard processes (i.e. from hadrons having pj > 
1 GeV/c) obtained from the scaled pQCD calculations (see 
later). Such a “perturbative” component accounts for a roughly 
constant ^7% factor of the total hadron multiplicity for all 
centralities. The good reproduction of the measured charged 
hadron integrated yields is an important result for our later use 
of dNch/dT\\x\=Q as an empirical measure of the initial entropy 
density in different Auh-Au centrality classes (see Section^}. 



Fig. 1. Charged hadron multiplicity at midrapidity (normalized by the 
number of participant nucleon pairs) as a function of centrality (given 
by the number of participants, Npart) measured in Au-l-Au at = 
200 GeV by PHENIX El (circles), STAR (42l (stars), PHOBOS El 
(squares) and BRAHMS 1441 (crosses), compared to our hydrodynam¬ 
ics calculations (dashed line), our scaled pQCD {pj > 1 GeV/c) p-l-p 
yields Ea (dashed-dotted line), and to the sum hydro-l-pQCD (solid 
line). 

For the initial conditions (Table E)> we choose ^o = 220 
GeV/fm^ (maximum energy density atb = 0 fm, corresponding 
to an average energy density over the total volume for head-on 
collisions of (eq) = 72 GeV/fm^) at a time Xq = 2R/y« 0.15 
fm/c equal to the transit time of the two Au nuclei at yTw = 200 
GeV. The choice of this relatively short value of Xq, - otherwise 
typically considered in other hydrodynamical studies of ther¬ 
mal photon production at RHIC II 011 41461 -, rather than the 
“standard” thermalization time of Ztherm = 0.6 fm/c II27I28I29I . 
is driven by our will to consistently take into account within our 
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Table 1. Summary of the thermodynamical parameters characterizing our hydrodynamical model evolution for central {b = Q fm) Au+Au 
collisions at = 200 GeV. Input parameters are the (maximum) initial energy density Eq (with corresponding ideal-gas entropy densities jq 
and temperature Tq) at time Xq, the baryochemical potential and the chemical and kinetic freeze-out temperatures T^kem and Tfg (or energy 
density Efo)- The energy densities at the end of the pure QGP (e'qqp), and at the beginning of the pure hadron gas phase (eJJJo) are also given for 
indication, as well as the average (over total volume) values of the initial energy density (eq), entropy density (sq), and temperature (7b). 


Xo 

(fm/c) 

eo ((eo)) 
(GeV/W) 

■so ((.so)) 
(fm-3) 

To ((7b)) 
(MeV) 

pinui 

^QGP 

(GeV/W) 

pmax 

^HRG 

(GeV/W) 

PB 

(MeV) 

T^chem 

(MeV) 

Tfo 

(MeV) 

^fo ^HRG 

(GeV/W) 

0.15 

220 (72) 

498 (190) 

590 (378) 

1.7 

0.35 

25. 

150 

120 

0.10 


space-time evolution the emission of photons from secondary 
“cascading” parton-parton collisions i47l4^ taking place in 
the thermalizing phase between prompt pQCD emission (at x ~ 
1/Pt ^0. 15 fm/c) and full equilibration (see Sect. 13.3t . Though 
it may be questionable to identify such photons from second- 
chance parton-parton collisions as genuine thermal y, it is clear 
that their spectrum reflects the momentum distribution of the 
partons during this thermalizing phase®. Additionally, recent 
theoretical works l49l^ do seem to support the application of 
hydrodynamics equations in such “pre-thermalization” condi¬ 
tions. Our consequent space-time evolution leads to a value of 
the energy density of e « 30 GeV/fm^ at ithenn = 0.6 fm/c, in 
perfect agreement with other 2 Dh- 1 hydrodynamic calculations 
which do not invoke azimuthal symmetry f!7rsi as well as 
more numerically involved 3 Dh- 1 approaches |2^. Thus, our 
calculations reproduce the final hadron spectra as well, at least, 
as those other works do. As a matter of fact, by using Tq = 0. 15 
fm/c (rather than 0.6 fm/c), the system has a few more tenths of 
fm/c to develop some extra transverse collective flow and there 
is no need to consider in our initial conditions a supplemen¬ 
tal input radial flow velocity parameter, as done in other 
works in order to reproduce the hadron spectra. 


2.3 Comparison to hadron data 

Figure |2] shows the pion, kaon, and proton^ transverse spec¬ 
tra measured by PHENIX Ea, STAR iEED, PHOBOS El 
and BRAHMS 1^ in central (0-10% corresponding to {b) = 
2.3 fm) and peripheral (60-70% corresponding to {b) = 11.9 
fm) Auh-Au collisions at = 200 GeV, compared to our hy¬ 
drodynamical predictions (dashed lines) and to properly scaled 
pH-p NLO pQCD expectations |45| (dotted lines). At low trans¬ 
verse momentum, the agreement data-hydro is excellent start¬ 
ing from the very low pj PHOBOS data {pr <100 MeV) up 
to at least pr ~ 1.5 GeV/c. Above this value, contributions 
from perturbative processes (parton fragmentation products) 
start to dominate over bulk hydrodynamic production. Indeed, 
particles with transverse momenta pr ^ 2 GeV/c are mostly 
produced in primary parton-parton collisions at times of or¬ 
der X ~ 1/pr ^ 0.15 fm/c (i.e. during the interpenetration of 

® Note also that it is precisely those secondary paitonic interactions 
that are actually driving the system towards (local) thermal equilib¬ 
rium. 

^ For a suitable comparison to the (feed-down corrected) 
PHENIX (52), PHOBOS (55l and BRAHMS E) yields, the STAR 
proton spectra 1531 have been appropriately corrected for a ~40% 
(p 7 -independent) contribution from weak decays 1541 . 


the colliding nuclei and before any sensible time estimate for 
equilibration), and as such, they are not in thermal equilibrium 
with the bulk particle production. Therefore, one does not ex¬ 
pect hydrodynamics to reproduce the spectral shapes beyond 
Pt ^ 2 GeV/c. The dotted lines of Fig. |2] show NLO pre¬ 
dictions for 7t, K and p production in pH-p collisions at y/s = 
200 GeV l45l scaled by the number of point-like collisions 
(Ncoii Taa) times an empirical quenching factor, Raa = 0.2 
(0.7) for 0-10% central (60-70% peripheral) Auh-Au, to account 
for the observed constant suppression factor of hadron yields at 
high Pt 1561571 (such a suppression is not actually observed in 
the p,p spectra at intermediate pT ^ 3 - 5 GeV/c, see discus¬ 
sion below). 

Fig. a shows more clearly (in linear rather than log scale 
as the previous figure) the relative agreement between the ex¬ 
perimental hadron transverse spectra and the hydrodynamical 
plus (quenched) pQCD yields presented in this work. The data- 
over-theory ratio plotted in the figure is obtained by taking the 
quotient of the pion, kaon and proton data measured in cen¬ 
tral Auh-Au reactions (shown in the left plot of Fig. |2j over the 
corresponding sum of hydrodynamical plus perturbative results 
(solid lines in Fig.|2j. In the low pT range dominated by hydro¬ 
dynamical production, there exist some local pj-dependent de¬ 
viations between the measurements and the calculations. How¬ 
ever, the same is true within the independent data sets them¬ 
selves and, thus, those differences are indicative of the amount 
of systematic uncertainties associated with the different mea¬ 
surements. High Pt hadro-production, dominated by perturba¬ 
tive processes, agrees also well within the ~20% errors as¬ 
sociated with the standard scale uncertainties for pQCD cal¬ 
culations at this center-of-mass energy. It is, thus, clear from 
Figs. m and 13 that identified particle production at y = 0 in 
nucleus-nucleus collisions at RHIC can be fully described in 
their whole pT range and for all centralities by a combination 
of hydrodynamical (thermalH-collective boosted) emission plus 
(quenched) prompt perturbative production. An exception to 
this rule are, however, the (anti)protons 1581 . Although due to 
their higher masses, they get an extra push from the hydro- 
dynamic flow up to Pr ^ 3 GeV/c, for even higher transverse 
momenta the combination of hydro plus (quenched) pQCD still 
clearly undershoots the experimental proton spectra. This ob¬ 
servation has lent support to the existence of an additional mech¬ 
anism for baryon production at intermediate pr values (pr « 
3-5 GeV/c) based on quark recombination l59l . This mecha¬ 
nism will not, however, be further considered in this paper since 
it has no practical implication for photon production and/or for 
the overall hydrodynamical evolution of the reaction. The over- 
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Au+Au (peripheral) calculations: 

. Hydro r [60-70%] 

. Hydro K [60-70%] xIO^ 

. Hydrop[60-70%]x10' 

. [NLOpQCD]xT^[60-70%],ii° 

. [NLOpQCD]xT„[60-70%],K'x10'' 

. [NLOpQCD]xT„[60-70%),px10'' 


AufAu (peripheral) data: 

• 

PHENIX JT [60-70%] 

• 

PHENIX Jt° [60-70%] 

0 

STARr [60-70%] 

■ 

PHENIXK* [60-70%]x10" 

☆ 

STARK* [60-70%] X 10’ 


stark: [60-80%] X 10’ 

k 

PHENIX p [60-70%] X 10^ 

☆ 

STAR p [60-70%] X 10"^ 



4 5 

Pt (GeV/c) 


4 5 

Pt (GeV/c) 


Fig. 2. Transverse momentum spectra for and protons measured in the range pj = 0 - 5.5 GeV/c by PHENIX 1^ . STAR {K° are 

preliminary) I53I54I . PHOBOS 1551 and BRAHMS 1311 in central (0-10% centrality, left) and peripheral (60-70%, right) Au-l-Au collisions at 
= 200 GeV, compared to our hydrodynamics calculations (dashed lines), to the scaled pQCD p-l-p rates 1451 (dotted lines), and to the sum 
hydro-l-pQCD (solid lines). 


all good theoretical reproduction of the differential %,K^p ex¬ 
perimental spectra for all centralities is obviously consistent 
with the previous observation that our calculated total inte¬ 
grated hadron multiplicities agree very well with the exper¬ 
imental data measured at mid-rapidity by the four different 
RHIC experiments (Fig.[3. 


3 Direct photon production 

As in the case of hadron production, the total direct photon 
spectrum in a given Au-i-Au collision at impact parameter b 
is obtained by adding the primary production from perturba¬ 
tive parton-parton scatterings to the thermal emission rates in¬ 
tegrated over the whole space-time volume of the produced 
fireball. Three sources of direct photons are considered cor¬ 
responding to each one of the phases of the reaction: prompt 
production, partonic gas emission, and hadronic gas radiation. 


3.1 Prompt photons 

For the prompt y production we use the NLO pQCD pre¬ 
dictions of W. Vogelsang i60l scaled by the corresponding 
Glauber nuclear overlap function at b, TAAib), as expected 
for hard processes in A-i-A collisions unaffected by final-state 
effects (as empirically confirmed for photon production in 
Au-i-Au II23I ). This pQCD photon spectrum is obtained with 
CTEQ6M 1611 parton distribution function (PDF), GRV 1621 
parametrization of the q,g ^ J fragmentation function (FF), 


and renormalization-factorization scales set equal to the trans¬ 
verse momentum of the photon (p = pr). Such NLO calcula¬ 
tions provide an excellent reproduction of the inclusive direct 
Y 1^ and large-p 7 7t® ED spectra measured by PHENIX in 
p-i-p collisions at = 200 GeV without any additional pa¬ 
rameter (in particular, at variance with results at lower ener¬ 
gies tm, no primordial kr is needed to describe the data). We 
do not consider any modification of the prompt photon yields 
in Au-i-Au collisions due to partially counteracting initial-state 
(IS) effects such as: (i) nuclear modifications (“shadowing”) of 
the Au PDE (< 20%, in the relevant (x,Q^) kinematical range 
considered here I46I21I63I ). and (ii) extra nuclear kj broaden¬ 
ing (Cronin enhancement) as described e.g. in 1201 . Both IS ef¬ 
fects are small and/or approximately cancel each other at mid¬ 
rapidity at RHIC as evidenced experimentally by the barely 
modified nuclear modification factor, RdAu for Y and TtP 

measured in d-i-Au collisions at = 200 GeV 1641 . Like¬ 
wise, we do not take into account any possible final-state (ES) 
photon suppression due to energy loss of the jet-fragmentation 
(aka. “anomalous”) component of the prompt photon cross- 
section I20I63I65I67I . which, if effectively present (see (511 
and discussion in Sect. 1.3.3> . can be in principle experimentally 
determined by detailed measurements of the isolated and non¬ 
isolated direct photon baseline spectra in p-i-p collisions at ^Js 
= 200 GeV |62l. 


3.2 Thermal photon rates 

Eor the QGP phase we use the most recent full leading order (in 
ttem and (Xs couplings) emission rates from Arnold et al. 1691 . 
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Au-fAu (central) data: 

• PHENIX 7C- [0-5%] 

• PHENIX Tco [0-10%] 

O STAR 71^5-10%] 

^ PHOBOS7t[0-15%] 
0 BRAHMS n- [0-5%] 
■ PHENIX K" [0-5%] 
☆ STAR [5-10%] 

■ir STARK°[0-5%] 

^ PHOBOSK[0-15%] 
A BRAHMS K' [0-5%] 
A PHENIX p [0-5%] 
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Fig. 3. Ratio of and protons yields measured in the range pj- = 0 - 5.5 GeV/c by PHENIX 1521 . STAR (note that are prelimi¬ 

nary) I53I54I . PHOBOS 1551 and BRAHMS l3ll in 0-10% most central Au-l-Au collisions at yTW = 200 GeV, over the sum hydro {{b) = 2.3 
fm) plus (quenched) pQCD. Theoretical calculations above pj ~2 GeV/c have an overall ±20% uncertainty (not shown) dominated by pQCD 
scale uncertainties. 


These calculations include hard thermal loop diagrams to all 
orders and Landau-Migdal-Pomeranchuk (LPM) medium in¬ 
terference effects. The parametrization given in 1169 1 assumes 
zero net baryon density (i.e. null quark chemical potential, 

= 0), and chemical together with thermal equilibrium. Cor¬ 
rections of the QGP photon rates due to net quark densities 
are o\py{%Tf] |20| i.e. marginal at RHIC energies where 
the baryochemical potential is close to zero at midrapidity 
{pB — 3 pcj ^ 25 MeV) and neglected here. Similarly, although 
the early partonic phase is certainly not chemically equilibrated 
(the first instants of the reaction are strongly gluon-dominated) 
the two main effects from chemical non-equilibrium composi¬ 
tion of the QGP: reduction of quark number and increase of the 
temperature, nearly cancel in the photon spectrum fTlTTl and 
have not been considered either. For the HRG phase, we use the 
latest improved parametrization from Turbide et al. EH which 
includes hadronic emission processes not accounted for be¬ 
fore in the literature. In all calculations, we use a temperature- 
dependent parametrization of the strong coupling^, (Xs{T) = 
2.095/{^ln(e/Aj^) ± [21n(e/Aj^)]} with Q = 2%T, 

obtained from recent lattice results l73l 


3.3 Extra photon contributions 

Apart from the aforementioned (prompt and thermal) photon 
production mechanisms, S. Bass et al. 14714^ have recently 
evaluated within the Parton Cascade Model (PCM), the con¬ 
tribution to the total Au±Au photon spectrum from secondary 
(cascading) parton-parton collisions taking place before the at¬ 
tainment of thermalization (i.e. between the transit time of the 
two nuclei, x « 0.15 fm/c, and the standard itherm = 0.6 fm/c 
considered at RHIC). Since such cascading light emission is 

® According to this parametrization, Ous{T) = 0.3 - 0.6 in the range 
of temperatures of interest here {T ~ 600 - 150 MeV). 


due to second-chance partonic collisions which are, simultane¬ 
ously, driving the system towards equilibrium, we consider not 
only “valid” (see the discussion of refs. |49|50|) but more self- 
consistent within our framework to account for this contribu¬ 
tion with our hydrodynamic al evolution alone. We achieve this 
by starting hydrodynamics (whose photon rates also include 
the expected LPM reduction of the secondary rates 1481 ) at Xo 
= 0.15 fm/c. By doing that, at the same time that we account 
for this second-chance emission, our initial plasma temperature 
and associated thermal photon production can be considered to 
be at their maximum values for RHIC energies. 


Likewise, we do not consider the conjectured extra y emis¬ 
sion due to the passage of quark jets (Compton-scattering and 
annihilating) through the dense medium 174I75I68I since such 
contribution is likely partially compensated by: (i) the con¬ 
current non-Abelian energy loss of the parent quarks going 
through the system 1661 . plus (ii) a possible photon suppres¬ 
sion due to energy loss of the “anomalous” component of the 
prompt photon cross-section 1201631671651 . As a matter of fact, 
some approximate cancellation of all those effects must ex¬ 
ist since the experimental Au±Au photon spectra above pr « 

4 GeV/c turn out to be well reproduced by primary (pQCD) 
hard processes alone for all centralities, as can be seen in the 
comparison of pQCD NLO predictions with PHENIX data l23l 
(Fig.S . The apparent agreement between the experimental spec¬ 
tra above pj k, A GeV/c and the NLO calculations does not 
seem to leave much room for extra radiation contributions. A 
definite conclusion on the existence or not of FS effects on 
photon production will require in any case precision y data in 
Au±Au, d±Au and p±p collisions. The more critical issue of 
the role of the jet bremsstrahlung component needs to be es¬ 
timated, for example, via measurements of isolated and non¬ 
isolated direct photon baseline p±p spectra as discussed in Ell. 
Additional IS effects not considered so far due, for example, to 
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isospin corrections® will require a careful analysis and compar¬ 
ison of Auh-Au to reference dn-Au photon cross-sections too. 


3.4 Total direct photon spectra 

Figure shows our computed total direct photon spectra for 
central (left) and peripheral (right) Auh-Au collisions at 
= 200 GeV, with the pQCD, QGP, and HRG components dif¬ 
ferentiated^* **. In central reactions, thermal photon production 
(mainly of QGP origin) outshines the prompt pQCD emission 
below pt ^ 3 GeV/c. Within pj- ~ 1-4 GeV/c, thermal pho¬ 
tons account for roughly 90% - 50% of the total photon yield in 
central Auh-Au, as can be better seen in the ratio total-y/pQCD- 
Y shown in Fig.|3 Photon production in peripheral collisions is, 
however, clearly dominated by the primary parton-parton radi¬ 
ation. In both cases, hadronic gas emission prevails only for 
lower pt values. In Figure we also compare our computed 
spectra to the inclusive Auh-Au photon spectra published re¬ 
cently by the PHENIX collaboration 1^ . The total theoretical 
(pQCDH-hydro) differential cross-sections are in good agree¬ 
ment with the experimental yields, though for central reactions 
our calculations tend to “saturate” the upper limits of the data 
in the range below « 4 GeV/c where thermal photons dom¬ 
inate. New preliminary PHENIX Auh-Au direct-^ results lEH 
1781 are also systematically above (tough still consistent with) 
these published spectra in the range pj- « 1-4 GeV/c and, if 
confirmed, will bring our results to an even better agreement 
with the data. 


To better distinguish the relative amount of thermal radi¬ 
ation in the theoretical and experimental total direct photon 
spectra in central Auh-Au collisions, we present in Eigure|3the 
nuclear modification factor defined as the ratio of the to¬ 
tal over prompt (i.e. 7)r4-scaled pH-p pQCD predictions) photon 
yields: 


i^AAiPr) 


dNTAVdpT 

Taa ■ da],p^^°/dpT 


( 2 ) 


A value « 1 would indicate that all the photon yield can 
be accounted for by the prompt production alone. Of course, 
since our total direct-y result for central Auh-Au includes ther¬ 
mal emission from the QGP and HRG phases, we theoretically 
obtain « 10 - 1 in the pj- « 1-4 GeV/c region where the 
thermal component is significant (Pig.|5|l. In this very same pj 
range, although the available PHENIX results have still large 
uncertainties**, the central value of most of the data points is 
clearly consistent with the existence of a significant excess over 


® Direct photon cross-sections depend on the light quark electric 
charges and are thus disfavoured in a nucleus target less rich in up 
quarks than the standard proton reference Ga 
*** We split the mixed phase contrihution onto QGP and HRG com¬ 
ponents calculating the relative proportion of QGP (HRG) matter in 
it. 

** Technically, the PHENIX data points helow pj = A GeV/c have 
“lower errors that extend to zero”, i.e. a non-zero direct-y signal is 
indeed observed in the data hut the associated errors are larger than 
the signal itself (S). 



Fig. 5. Direct photon “nuclear modification factor”, (Eq.|^, ob¬ 
tained as the ratio of the total over the prompt y spectra for 0-10% 
most central Auh-Au reactions at = 200 GeV. The solid line is 
the ratio resulting from our hydroH-pQCD model. The points show the 
PHENIX data 1231 over the same NLO yields and the dashed-dotted 
curves indicate the theoretical uncertainty of the NLO calculations 
(see text). 


the NLO pQCD expectations. A note of caution is worth here, 
however, regarding the ^ 1 value observed for both the 
theoretical and experimental spectra below pj- « 4 GeV/c since 
it is not yet clear to what extent the NLO predictions, enter¬ 
ing in the denominator of Eq. 0 , are realistic in this thermal- 
photon “region of interest”. Indeed, in this comparatively low 
Pt range the theoretical prompt yields are dominated by the 
jet bremsstrahlung contribution El which is intrinsically non- 
perturbative (i.e. not computable) and determined solely from 
the parametrized parton-to-photon GRV 1621 EE which is rela¬ 
tively poorly known in this kinematic range. The standard scale 
uncertainties in the NLO pQCD calculations are ±20% above 
Pt K, A GeV/c but we have assigned a much more pessimistic 
^ 5 Q**% uncertainty to these calculations in the range pT ~ \ - 
4 GeV/c (dashed-dotted lines in Eig.|3. Precise measurements 
of the direct-y baseline spectrum in pH-p collisions at y/s = 200 
GeV above pj- = 1 GeV/c are mandatory before any definite 
conclusion can be drawn on the existence or not of a thermal 
excess from the Auh-Au experimental data. 

As a final cross-check of our computed hydrodynamical 
photon yields, we have compared them to previously published 
predictions for thermal photon production in Auh-Au collisions 
at top RHIC energy: D. K. Srivastava et al. cni (with initial 
conditions Tq « 0.2 fm/c and Tq « 450 - 660 MeV), Jan-e Alam 
et alP CD (70 = 0.5 fm/c and Tq = 300 MeV), E. D. Steffen 
and M. H. Thoma CD (To = 0.5 fm/c and 7b = 300 MeV), 
S. S. Rasanen et al. C3 (To = 0.17 fm/c and 7b = 580 MeV), 
N. Hammon etal. 1461 (Tq = 0.12 fm/c and To = 533 MeV), and 
Turbide et al.^^ ED (Xo = 0.33 fm/c and 7b = 370 MeV). Eor 

*^ Alam et al. have recently 1771 recomputed their hydrodynamical 
yields using higher initial temperatures (Tq = 400 MeV at Xq = 0.2 
fm/c) and getting a better agreement with the data. 

*^ Note that stricto senso Turbide’s spectra are not obtained with 
a pure hydrodynamical computation but using a simpler “dynamical 























8 


David d’Enterria, and Dmitri Peressounko: Probing the QCD EoS with thermal yin A+A collisions at RHIC 




Fig. 4. Photon spectra for central (0-10%, left) and peripheral (60-70%, right) Au+Au reactions at = 200 GeV as computed with our 
hydrodynamical model [with the contributions for the QGP and hadron resonance gas (HRG) given separately] compared to the expected NLO 
pQCD p+p yields for the prompt y 1601 (scaled by the corresponding nuclear overlap function), and to the experimental photon yields measured 
by the PHENIX collaboration 1231 . 



Fig. 6. Thermal photon predictions for central Au+Au reactions at 
= 200 GeV as computed with different hydrodynamical 1101 
I11I14I or “dynamical fireball” 1721 models, compared to (i) our hy¬ 
dro calculations (dashed curve), (ii) the expected perturbative y yields 
(7)iA-scaled NLO p+p calculations I60I L and (iii) the experimental to¬ 
tal direct photon spectrum measured by PHENIX 1231 . 


similar initial conditions, the computed total thermal yields in 
those works are compatible within a factor of ^2 with those 
presented here. Some of those predictions are shown in Fig¬ 
ure confronted to our calculations. Our yields are, in gen¬ 
eral, above all other predictions since, as aforementioned, both 
our initial thermalization time and energy densities (tempera¬ 
tures) have the most “extreme” values possible consistent with 
the RHIC charged hadron multiplicities. They agree specially 
well with the hydrodynamical calculations of the Jyvaskyla 
group CH which have been computed with the same up-to- 
date QGP rates used here. Given the current (large) uncertain¬ 
ties of the available published data, all thermal photon predic¬ 
tions are consistent with the experimental results. However, as 
aforementioned, newer (preliminary) PHENIX direct-y* mea¬ 
surements have been reported very recently IPI781 and indi¬ 
cate a clear excess of direct photons over NLO pQCD for Au 
+Au at = 200 GeV in this pr range in excellent agree¬ 
ment with our thermal photon calculations. 

4 Thermal photons and the QCD 
equation-of-state 

In order to experimentally isolate the thermal photon spectrum 
one needs to subtract from the total direct y spectrum the non- 
equilibrated “background” of prompt photons. The prompt y 
contribution emitted in a given Au+Au centrality can be mea¬ 
sured separately in reference p+p (or d+Au) collisions at the 


fireball” model which assumes constant acceleration in longitudinal 
and transverse directions. 
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same scaled by the corresponding nuclear overlap function 
TAA{b), and subtracted from the total Au+Au y spectrum 1671 . 
The simpler expectation is that the remaining photon spectrum 
for a given impact parameter b 

dNZl'Jib) daZ 

-7- = -7- ^AA(b )--—, (3) 

dpT dpT dpT 

will be just that due to thermal emission from the partonic and 
hadronic phases of the reaction. Such a subtraction procedure 
can be effectively applied to all the y spectra measured in dif¬ 
ferent centralities as long as both the total Auh-Au and baseline 
pH-p photon spectra are experimentally measured with reason¬ 
able (<15%) point-to-point (systematical and statistical) uncer¬ 
tainties ia. The subtracted spectra 0 can be therefore sub¬ 
ject to scrutiny in terms of the thermodynamical properties of 
the radiating medium. 


4.1 Determination of the initial temperature 

Due to their weak electromagnetic interaction with the sur¬ 
rounding medium, photons produced in the reaction escape 
freely the interaction region immediately after their production. 
Thus, even when emitted from an equilibrated source, they are 
not reabsorbed by the medium and do not have a black-body 
spectrum at the source temperature. Nonetheless, since all the 
theoretical thermal y rates l6W7^ have a general functional 
dependence of the form^^ EydRy/d^p ■ exp{—Ey/T), 

one would expect the hnal spectrum to be locally exponential 
with an inverse slope parameter strongly correlated with the 
(local) temperature T of the radiating medium. Obviously, 
such a general assumption is complicated by several facts. 
On the one hand, the final thermal photon spectrum is a sum 
of exponentials with different temperatures resulting from 
emissions at different time-scales and/or from different regions 
of the fireball which has strong temperature gradients (the 
core being much hotter than the “periphery”). On the other 
hand, collective flow effects (stronger for increasingly central 
collisions) superimpose on top of the purely thermal emis¬ 
sion leading to an effectively larger inverse slope parameter 
{T,ff « v'(l+P)/(l”P)7') O- One of the main results of 
this paper is to show that, based upon a realistic hydrody- 
namical model, such effects do not destroy completely the 
correlation between the apparent photon temperature and the 
maximal temperature actually reached at the beginning of the 
collision process. We will show that such a correlation indeed 
exists and that the local inverse slope parameter obtained by 
htting to an exponential, at high enough pr, the thermal photon 
spectrum obtained via the expression 0, indeed provides a 
good proxy of the initial temperature of the system without 
much distortion due to collective flow (and other) effects. 

To determine to what extent the thermal slopes are indica¬ 
tive of the original temperature of the system, we have fitted 

The factor is just an overall normalization factor in this case 
(since its temporal variation is small compared to the short emission 
times) and does not significantly alter the exponential shape of the 
spectra. 


the thermal spectra obtained from our hydrodynamical calcula¬ 
tions in different Auh-Au centralities to an exponential distribu¬ 
tion in different pj ranges. Since, - according to our Glauber 
prescription for the impact-parameter dependence of the hy¬ 
drodynamical initial conditions -, different centralities result 
in different initial energy densities, we can in this way explore 
the dependence of the apparent thermal photon temperature on 
the maximal initial temperatures Tq (at the core) of the sys¬ 
tem. The upper plot of Figure0shows the obtained local slope 
parameter, T^jj, as a function of the initial To for our default 
QGPh-HRG hydrodynamical evolution (Table 0. We find that 
although all the aforementioned effects smear the correlation 
between the apparent and original temperatures, they do not 
destroy it completely. The photon slopes are indeed approxi¬ 
mately proportional to the initial temperature of the medium, 
Tq. There is also an obvious anti-correlation between the pj 
of the radiated photons and their time of emission. At high 
enough pj the hardest photons issuing from the hottest zone 
of the system swamp completely any other softer contributions 
emitted either at later stages and/or from outside the core region 
of the fireball. Thus, the higher the pj range, the closer is Tgjf 
to the original Tq at the center of the system. According to our 
calculations, empirical thermal slopes measured above pj thA 
GeV/c in central Auh-Au collisions are above ~400 MeV i.e. 
only ~30% lower than the “true” maximal (local) temperature 
of the quark-gluon phase. On the other hand, local y slopes 
in the range below pj- « 1 GeV/c have almost constant value 
Teff ^ 200 MeV (numerically close to 7)r,f) for all centralities 
and are almost insensitive to the initial temperature of the hy¬ 
drodynamical system but mainly specified by the exponential 
prefactors in the hadronic emission rates, plus collective boost 
effects. 

To assess the dependence of the thermal photon spectra on 
the underlying EoS, we have rerun our hydro evolution with 
just the EoS of a hadron resonance gas. We choose now as ini¬ 
tial conditions: eo = 30 GeV/fm^ at Tq = 0.6 fm/c, which can 
still reasonably describe the experimental hadron spectra. Ob¬ 
viously, any description in terms of hadronic degrees of free¬ 
dom at such high energy densities is unrealistic but we are in¬ 
terested in assessing the effect on the thermal photon slopes 
of a non ideal-gas EoS as e.g. that of a HRG-like system with 
a large number of heavy resonances (or more generally, of any 
EoS with exponentially rising number of mass states). The pho¬ 
ton slopes for the pure HRG gas EoS (Eig.0 bottom) are lower 
(r^J^ « 220 MeV) than in the default QGPh-HRG evolution, 
not only because the input HRG Eq is smaller (the evolution 
starts at a later Tq) but, specially because for the same ini¬ 
tial Eo the effective number of degrees of freedom in a system 
with a HRG EoS is higher than that in a QGP^^ and therefore 
the initial temperatures are lower. A second difference is that, 
for all pt ranges, we find almost the same exact correlation 
between the local y slope and To indicating a single underly¬ 
ing (hadronic) radiation mechanism dominating the transverse 
spectra at all pr- 


Note that g{T) e/T^ increases exponentially with T for a HRG- 
like EoS, and at high enough temperatures will clearly overshoot the 
QGP constant number of degrees of freedom. 
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Fig. 7. Local photon slope parameters Tgff (obtained from exponential 
fits of the thermal photon spectrum in different pj ranges) plotted 
versus the initial (maximum) temperature Tq of the fireball produced 
at different centralities in Au+Au collisions at ^/s = 200 GeV. Upper 
plot - hydrodynamical calculations with QGP+HRG EoS (Table 0, 
bottom - HRG EoS (with initial conditions: Eq = 30 GeV/fm^ at Tq = 
0.6 fm/c). 


Two overall conclusions can be obtained from the study of 
the hydrodynamical photon slopes. First, the observation in the 
data, via Eq. (|3, of a thermal photon excess above pr « 2.5 
GeV/c with exponential slope Tgjj > 250 MeV is an unequivo¬ 
cal proof of the formation of a system with maximum temper¬ 
atures above Tcrit since no realistic collective flow mechanism 
can generate such a strong boost of the photon slopes, while 
simultaneously reproducing the hadron spectra. Secondly, pro¬ 
nounced pt dependences of the local thermal slopes seem to 
be characteristic of space-time evolutions of the reaction that 
include an ideal-gas QGP radiating phase. 


is implicit in our perfect fluid hydrodynamical equations with 
zero viscosity) one can estimate the initial entropy density s at 
the time of photon emission from the total final particle multi¬ 
plicity dN/dy measured in the reaction. Varying the centrality 
of the collision, one can then explore the form of the depen¬ 
dence s = s{T) at the first instants of the reaction, extract the 
underlying equation of state of the radiating system and trace 
any signal of a possible phase transition. Indeed, the two most 
clear evidences of QGP formation from QCD calculations on 
the lattice are: (i) the sharp rise of e{T)/T^, or equivalently of 
s{T)/T^, at temperatures around Tcrit, and (ii) the flattening of 
the same curve above Tcrit- The sharp jump is of course due to 
the sudden release of a large number of (partonic) degrees of 
freedom at Tcrit- The subsequent plateau is due to the full for¬ 
mation of a QGP with a fixed (constant) number of degrees of 
freedom. 

We propose here to use Tejf as a proxy for the initial tem¬ 
perature of the system, and directly study the evolution, versus 
Tgff, of the effective number of degrees of freedom defined as^® 

^ ^ ^ ( 4 ) 

which coincides with the degeneracy of a weakly interacting 
gas of massless particles. [In a similar avenue, B. Muller and 
K. Rajagopal (Si have recently proposed a method to esti¬ 
mate the number of thermodynamic degrees of freedom via 
geff /^, where s is also determined from the final hadron 
multiplicities]. The dashed line in Fig. El (top) shows the evo¬ 
lution of the true number of degrees of freedom giiydro{so,To) 
computed via Eq. 0 , as a function of the (maximal) temper¬ 
atures and entropies directly obtained from the initial condi¬ 
tions of our hydrodynamical model in different Auh-Au cen¬ 
tralities'^. The first thing worth to note is that g{s,T) remains 
constant at the expected degeneracy ghydro = 42.25 of an ideal 
gas of Nf = 2.5 quarks and gluons for basically all the max¬ 
imum temperatures accessible in the different centralities of 
Auh-Au at yTw = 200 GeV. This indicates that at top RHIC 
energies and for most of the impact parameters. To is (well) 
above Tcrit and the hottest parts of the initial fireball are in the 
QGP phase. The expected drop in ghydro related to the transi¬ 
tion to the hadronic phase is only seen, if at all, for the very 
most peripheral reactions (with To ~ Tc)- Thus, direct evidence 
of the QGP-HRG phase change itself via the study of the cen¬ 
trality dependence of any experimentally accessible observable 
would only be potentially feasible at RHIC in Auh-Au reactions 
at lower center-of-mass energies 1801 . 


4.2 Determination of the QCD Equation of State (EoS) 

As we demonstrated in the previous section, Tgff is approxi¬ 
mately proportional to the maximum temperature reached in a 
nucleus-nucleus reaction. One can go one step further beyond 
the mere analysis of the thermal photon slopes and try to get a 
direct handle on the equation of state of the radiating medium 
by looking at the correlation of Tejf with experimental observ¬ 
ables related to the initial energy or entropy densities of the 
system. For example, assuming an isentropic expansion (which 


As aforementioned, we can empirically trace the QCD EoS 
shown in Fig. |3 (and eventually determine the temperature- 
evolution of the thermodynamic degrees of freedom of the pro¬ 
duced medium) using the estimate of the initial temperature 
given by the thermal photon slopes, and a second observ¬ 
able closely related to the initial entropy of the system such as 

Units are in GeV and fm. ^(4) = 71^/90, where ^(n) is the Riemann 
zeta function. 

In the most peripheral reactions, the bag entropy has been sub¬ 
tracted to make more apparent the drop near fi- 
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Fig. 8. Effective initial number of degrees of freedom obtained from 
our hydrodynamical calculations with a QGP+HRG EoS (upper plot), 
and with a pure HRG EoS (bottom), plotted as a function of the tem¬ 
perature (Tq) or thermal photon slope (Tgff) in different Au-l-Au cen¬ 
trality classes at y/s^ = 200 GeV. The number of degrees of free¬ 
dom are computed respectively: (i) From our initial thermodynamical 
conditions (^q, To) via Eq. (dashed line), (ii) from the obtained 
charged hadron multiplicity dNch/dr\ and the true initial temperature 
Tq via Eq. (dotted-dashed line); and (iii) from dNch/dr[ and the 
thermal photon slopes Tgff measured in different pj ranges via Eq. 
(solid lines). For illustrative purposes, the open squares indicate the 
approximate position of the different Au-l-Au centrality classes (in 
10% percentiles) for the values of ggg obtained using the thermal pho¬ 
ton slopes measured above pr = 4 GeV/c. 


the final-state hadron multiplicity, dN/dy. Although one could 
have also considered to obtain ggjp via e/T^ {dEj/dy)/T^^, 
using the transverse energy per unit rapidity dEj /dy measured 
in different Auh-Au centralities ED, we prefer to use the ex¬ 
pression 0 which contains the entropy-, rather than the energy- 
, density for two reasons: 

(i) the experimentally accessible values of r/A/r/y remain con¬ 
stant in an isentropic expansion (i.e. dN/dy sq) whereas, 
due to longitudinal work, the measured final dEj/dy pro¬ 
vides only a lower limit on the initial e {dEj/dy < Eq); and 

(ii) geff is less sensitive to experimental uncertainties 

associated to the measurement of Tgg than ggjp 

Again, in the absence of dissipative effects, the space-time 
evolution of the produced system in a nucleus-nucleus reac¬ 
tion is isentropic and the entropy density (per unit rapidity) 


at the thermalization time To can be directly connected (via 
i « 4p 1811 ) to the final charged hadron pseudo-rapidity den¬ 
sity 

■ dN ^ 7.2 dNch 

dV (Ax) - To ^fr| 

where we have written the volume of the system, dV = 
(Ax)xoc/ri, as the product of the (purely geometrical) average 
transverse overlap area for each centrality times the starting 
proper time of our hydro evolution (xq =0.15 fm/c), and where 
dNch/dr\ is the charged hadron multiplicity customarily mea¬ 
sured experimentally at mid-rapidity^^. By combining, Eqs. 0 
and 0 , we obtain an estimate for the number of degrees of 
freedom of the system produced in a given Ah-A collision at 
impact parameter b: 


geff 


/ dNchjb) 

V 



150 (Sc)3 dNch{b) 


(6) 

which can be entirely determined with two experimental 
observables: dNch/dr\ and Tgff. 


Let us first assess to what extent the ansatz 0 is affected 
by the assumption that Eq. © indeed provides a good exper¬ 
imental measure of the initial entropy density s. The dotted- 
dashed curve in Eig. 0 has been obtained via Eq. 0 using 
the {dNch/dr\)/ {A\/) values obtained from our hydrodynam¬ 
ical model, and the true (input) initial temperature of the sys¬ 
tem To, and thus it is only sensitive to the way we estimate 
the entropy density. The resulting curve is a factor of ^3 be¬ 
low the expected “true” ghydm curve, i.e. g^g {dNch/dr[,To) ~ 
3 ■ g hydro i^o, To), indicating that Eq. 0 underestimates by the 
same amount the maximal entropy of the original medium. This 
is so because our estimate {dNch/dr[)/ {A±) specifies the en¬ 
tropy density averaged over the whole Glauber transverse area 
(Ax), whereas the maximal entropy area in the core of the sys¬ 
tem (from where the hardest thermal photons are emitted) is 
~3 times smaller. Although one could think of a method to 
correct for this difference, this would introduce an extra model- 
dependence that we want to avoid at this point. We prefer to 
maintain the simple (geometrical overlap) expression of the 
transverse area (Ax(^)) in Eq. 0 , and exploit the fact that, 
although such an equation does not provide the true absolute 
number of degrees of freedom, it does provide a very reliable 
indication of the dependence of ggjp on the temperature of the 
system and, therefore, of the exact/orm of the underlying EoS. 


Einally, let us consider the last case where we use Eq. 0 
with the values of dNch/dr\ and Teff that can be actually ex¬ 
perimentally measured. The different solid curves in the upper 
plot of Eig.l^show the effective degeneracy geg, computed us¬ 
ing Eq. 0 and the local photon slopes Teff measured in dif¬ 
ferent pt ranges for our default QGPh-HRG evolution. As one 
could expect from Eig. 0 the best reproduction of the shape 

This formula uses Ntot/Nch = 3/2, and the Jacobian |iir|/dy| = 
E/py=y 1 . 2 . 

Note again that both the photon slopes and the charged hadron 
multiplicities are proxies of the thermodynamical conditions of the 
system at the same time Xq. 
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of the underlying EoS is obtained with the effective tempera¬ 
tures measured in higher pj bins. For those Tejf, the computed 
ge^’s show a relatively constant value in a wide range of cen¬ 
tralities as expected for a weakly interacting QGP. Deviations 
from this ideal-gas plateau appear for more central collisions, 
due to an increasing difference between the (high) initial tem¬ 
peratures, 7b, and the apparent temperature given by the photon 
slopes (Fig.E). Such deviations do not spoil, however, the use¬ 
fulness of our estimate since, a non-QGP EoS would result in a 
considerably different dependence of geff on the reaction cen¬ 
trality. Indeed, the different curves in the bottom plot of Fig. |3 
obtained with a pure hadron resonance gas EoS clearly indi¬ 
cate^*' that a HRG EoS, or in general any EoS with exponen¬ 
tially increasing number of mass states, would bring about a 
much more dramatic rise of gejj with Tejf- 

In summary, the estimate (|6} indeed provides a direct ex¬ 
perimental handle on the form of the EoS of the strongly in¬ 
teracting medium produced in the first instants of high-energy 
nuclear collisions. More quantitative conclusions on the possi¬ 
bility to extract the exact shape of the underlying EoS and/or 
the absolute number of degrees of freedom of the produced 
medium require more detailed theoretical studies (e.g. with vary¬ 
ing lattice-inspired EoS’s jSOl and/or using more numerically 
involved 3 Dh- 1 hydrodynamical approaches). In any case, we 
are confident that by experimentally measuring the thermal pho¬ 
ton slopes in different Auh-Au centralities and correlating them 
with the associated charged hadron multiplicities as in Eq. 
one can approximately observe the expected “plateau” in the 
number of degrees of freedom indicative of QGP formation 
above a critical value of T. 


5 Conclusions 

We have studied thermal photon production in Auh-Au reac¬ 
tions at = 200 GeV using a Bjorken hydrodynamic model 
with longitudinal boost invariance. We choose the initial con¬ 
ditions of the hydrodynamical evolution so as to efficiently re¬ 
produce the observed particle multiplicity in central Auh-Au 
collisions at RHIC and use a simple Glauber prescription to 
obtain the corresponding initial conditions for all other central¬ 
ities. With such a model we can perfectly reproduce the iden¬ 
tified soft pion, kaon and proton pj-differential spectra mea¬ 
sured at RHIC. Complementing our model with the most up- 
to-date parametrizations of the QGP and HRG thermal photon 
emission rates plus a NLO pQCD calculation of the prompt y 
contribution, we obtain direct photon spectra which are in very 
good agreement with the Auh-Au direct photon (upper limit) 
yields measured by the PHENIX experiment. In central colli¬ 
sions, a thermal photon signal should be identifiable as a factor 
of ~8 - 1 excess over the pQCD y component within pr ~ I - 
4 GeV/c, whereas pure prompt emission clearly dominates the 
photon spectra at all pj in peripheral reactions. The local in¬ 
verse slope parameter of the thermal photon spectrum is found 
to be directly correlated to the maximum temperature attained 

Accidentally, > ghydm in the case of a HRG EoS, because 
the underestimation of the apparent temperature (raised to the cube) 
compensates for the aforementioned area averaging of the entropy. 


in the course of the collision. The experimental measurement 
of local thermal photon slopes above pj « 2.5 GeV/c, with 
values Tgff > 250 MeV and with pronounced pj dependences 
can only be reproduced by space-time evolutions of the reac¬ 
tion that include a QGP phase. 

Finally, we have proposed and tested within our framework, 
an empirical method to determine the effective thermodynam¬ 
ical number of degrees of freedom of the produced medium, 
s{T)/T^, by correlating the thermal photon slopes 
with the final-state charged hadron multiplicity measured in 
different centrality classes. We found that one can clearly dis¬ 
tinguish between the equation of state of a weakly interacting 
quark-gluon plasma and that of a system with rapidly rising 
number of mass states with T. Stronger quantitative conclu¬ 
sions on the exact shape of the underlying EoS and/or the abso¬ 
lute number of degrees of freedom of the produced medium re¬ 
quire more detailed theoretical studies as well as high precision 
photon data in Auh-Au and baseline pH-p, dn-Au collisions. In 
any case, the requirement for hydrodynamical models of con¬ 
currently describing the experimental bulk hadron and thermal 
photon spectra for different Auh-Au centralities at = 200 
GeV, imposes very strict constraints on the form of the equa¬ 
tion of state of the underlying expanding QCD matter produced 
in these reactions. 
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